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Abstract 

In this paper, an integrated flight/propulsion control system design is presented for 
the piloted longitudinal landing task with a modern, statically unstable, fighter aircraft. A 
centralized compensator based on the Linear Quadratic Gaussian/Loop Transfer Recovery 
methodology is first obtained to satisfy the feedback loop performance and robustness 
specifications. This high-order centralized compensator is then partitioned into airframe 
and engine "sub-controllers" based on modal controllability /observability for the 
compensator modes. The order of the sub— controllers is then reduced using 
internally-balanced realization techniques and the sub-controllers are simplified by 
neglecting the insignificant feedbacks. These sub-controllers have the advantage that they 
can be implemented as separate controllers on the airframe and the engine while still 
retaining the important performance and stability characteristics of the full-order 
centralized compensator. Command prefilters are then designed for the closed— loop system 
with the simplified sub— controllers to obtain the desired system response to airframe and 
engine command inputs, and the overall system performance evaluation results are 
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presented. 


Introduction 

The desire to enlarge the flight envelope of tactical aircraft and provide 
new/enhanced maneuver capabilities has led to the use of forces and moments produced by 
the propulsion system to augment the flight control function. The coupling between the 
propulsion system dynamics and the airframe dynamics, which results from such an 
application of the propulsion system forces and moments, is significant enough that the 
traditional approach of designing the propulsion control system and the flight control 
system separately and then putting them together in an ad— hoc manner is no longer 
adequate. An integrated approach to flight/propulsion control system design is then 
required in order to obtain an overall system that will ensure optimum performance with 
minimal pilot workload. 

In the early 1980s, the Air Force initiated a study called the Design Methods for 
Integrated Control Systems (DMICS) with the objective of developing methodologies for 
the design of Integrated Flight/Propulsion Control (IFPC) laws for an advanced tactical 
aircraft. Two very different approaches to IFPC design came about as a result of this 
study. These are : 

(i) A Linear Quadratic Gaussian — Loop Transfer Recovery (LQG/LTR) based 
centralized approach which consists of designing a "global" integrated compensator 
considering the fully integrated system as one high-order system W. 

(ii) A decentralized, hierarchical approach which consists of partitioning the integrated 
system into loosely coupled sub-systems and then designing separate controllers for 
the sub— systems, using Linear Quadratic Regulator (LQR) based explicit 
model-following type of approach, such that some high-level performance criterion 
are met^. 

More recently, a joint NASA Lewis and NASA Ames research program, called the Short 


2 



Take-Off and Vertical Landing (STOVL) Integrated Flight/Propulsion Controls 
program^, has been initiated with the objective of developing the enabling technologies 
that will lead to a demonstrator aircraft with Supersonic STOVL capability in the 1990s. 
One of the elements of this program is to develop IFPC design methodologies that will 
provide a viable alternative to the methodologies developed under the DMICS study while 
allowing for improved system performance and simplicity of control law synthesis and 
implementation. Before embarking on developing an alternative to DMICS, it is 
meaningful to conduct an extensive review of the two DMICS methodologies and develop a 
thorough understanding of the control law synthesis and evaluation procedures by way of 
an example study. Towards this goal, results are reported in this paper from an example 
application of the LQG/LTR— based DMICS approach to IFPC design for the piloted 
longitudinal landing task with a modern, statically unstable fighter aircraft. This study, 
however, takes the IFPC design process a step beyond that reported in Ref. [1], in that the 
high-order centralized compensator is partitioned into simplified lower-order 
sub— controllers which can be implemented separately on the airframe and the engine 
without any significant loss of overall system performance and robustness. 

In the following, the LQG/LTR— based DMICS control law design procedure is first 
briefly summarized. The vehicle model is then discussed, and the design specifications are 
stated. The centralized feedback compensator design is then presented and the closed-loop 
performance and robustness characteristics are discussed. The procedure for partitioning 
the centralized compensator into separate airframe and engine sub-controllers is discussed 
next and closed— loop evaluation results are presented for the system consisting of the 
simplified sub— controllers. Finally, prefilter design results are presented to provide 
command shaping such that the closed— loop system response to airframe and engine 
command inputs is as desired. 
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Control Law Design Methodology 

The salient features of the LQG/LTR— based DMICS methodology are briefly 
discussed in the following. This discussion is based on the description of the methodology in 
Ref. [4], 

The nonlinear simulation of the complete aircraft system is first used to obtain 
linear models at various design points along the aircraft flight envelope. An open-loop 
analysis of the plant, at each design point, is performed to assess the number of degrees of 
freedom and to define the candidate variables to be regulated and/or commanded. To 
determine the plant degrees of freedom, an input/output scaling (normalization) of the 
plant is performed. This scaling consists of normalizing the plant inputs by the maximum 
allowable deflection of the control effector and normalizing the plant outputs by the 
maximum value to be commanded or the maximum allowable deviation from trim. The 
singular values of the scaled plant transfer matrix plotted as a function of frequency then 
provide an indication as to whether adequate control authority is available to command the 
desired outputs by the available control inputs. At this stage, further dynamics are 
appended to the plant to provide desirable feedback— loop properties. An example of these 
augmentation dynamics is the addition of integrators in each control channel to ensure that 
the resultant system tracks step commands with zero steady— state error. 

The control design structure consists of feedback control with command shaping, as 
shown in the block diagram of Fig. 1. Design goals are used to generate requirements on 
feedback-loop properties and on command shaping prefilters. The feedback compensator, 
K(s), is designed first using a Linear Quadratic Gaussian with Loop Transfer Recovery 
(LQG/LTR)^’^ approach. The LQG/LTR methodology provides an integrated 
frequency-domain and state— space approach for design of multi— input multi— output 
(MIMO) control systems and directly addresses design issues such as stability robustness 
and trade-off between performance and allowable control power. The basic requirements 
that the closed— loop feedback system must satisfy are : ( Instability — bounded outputs, 
y(s), for all bounded disturbances, d(s), and bounded reference inputs, r(s), (2)Performance 
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- small errors, e(s), in the presence of r(s) and d(s), and (3)Robustness - stability and 
performance maintained in the presence of model uncertainties. These requirements can be 
transformed into specifications on the form of the loop shapes for the feedback system. The 
LQG/LTR procedure then consists of synthesizing the feedback compensator, K(s), by 
manipulating the LQG design parameters in such a way that the desired loop shapes are 
obtained. 

After the feedback design is complete, the command shaping, P(s), is designed with 
a multi-variable bandlimited inverse method^. This method is based on the theoretical 
asymptotic properties of the linear quadratic stochastic control problem and provides for 
setting the bandwidth of the prefiltering action by properly selecting the control weighting 
in the regulator problem. Such an approach allows the designer to obtain a prefilter, P(s), 
for desired response to command inputs without altering the feedback properties. 

The command shaping and feedback compensator designs can performed for the 
various design points covering the aircraft flight envelope and a practical controller for the 
full flight envelope obtained by gain scheduling the designs, adding limiting functions, etc.. 
In this study, however, the control design will be performed for only one fixed point along 
the approach and landing trajectory. 

A major criticism of a centralized IFPC design approach, such as the one discussed 
above, is that it results in one high order complicated compensator which is very difficult 
to implement. Such a compensator might also result in many feedback paths which are not 
physically realizable. Also, traditionally it is the responsibility of the engine 
designer/manufacturer to make sure that the engine will provide the desired performance 
when installed in the airframe. The engine manufacturer performs extensive "bench tests" 
with an active engine controller to assure an adequate design. The accountability issue, 
then, becomes a major political hurdle with a centralized airframe and engine compensator. 
To address these difficulties, the idea of partitioning the centralized compensator into 
separate airframe and engine sub-controllers, without any significant loss in overall system 
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performance, is introduced in this paper. The mathematics associated with the 
compensator partitioning and the example application are discussed in a later section. 


Vehicle Model 


The vehicle model consists of an integrated airframe and propulsion system 
state— space representation for a modern fighter aircraft powered by a two— spool turbofan 
engine and equipped with a 2D thrust— vectoring and reversing nozzle. The vehicle 
dynamics are linearized at an airspeed V Q = 120 Knots and flight path angle 7 Q = —3 deg. 
This flight condition is representative of the STOL (Short Take-off and Landing) 
approach— to-landing task. The vehicle model has the form 


x„ = A x„ -|- B,u ; y^ = C,x + D 11 
P PP PP ,J P PP PP 


where the state vector is 


Xp= [u, w, q, 6 , h, N2, N25, P6, T41B] 


with 


( 1 ) 


u = aircraft body axis forward velocity (ft /sec) 
w = aircraft body axis vertical velocity (ft/sec) 
q = aircraft pitch rate (rads/sec) 

0 = pitch angle (rads) 
h = altitude (ft) 

N2 = engine fan speed (rpm) 

N25 = core compressor speed (rpm) 

P6 = engine mixing plane pressure (psia) 

T41B = engine high pressure turbine blade temp. (°R), 
and the control input vector is 

Up= [<$pp, ^p'p) WF, A78, A8, <!tpy] 


with 


<!>prp = leading edge flap deflection angle(deg) 

= trailing edge flap deflection angle (deg) 

WF = engine main burner fuel flow rate (#/hr) 

o 

A78 = thrust reverser port area (in ) 

2 

A8 = main nozzle throat area (in ) 
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<5rpy = nozzle thrust vectoring angle (deg). 

The vector yp represents the outputs to be controlled and will be discussed later in this 
paper. The vehicle system matrix Ap and the control effectiveness matrix Bp are listed in 
the Appendix. 

The open— loop vehicle eigenvalues are : 

Aj= 0.07, A 2 2 = -0.09±j0.23, A^= 1.06, A^= — 1.47 — Airframe modes 

Ag= —1.40, \j= —3.57, Ag= -6.96, A^= —89.28 — Propulsion modes. 

Note that the airframe is statically unstable with a highly unstable pitch mode. The 

open— loop plant Ap matrix is full, i.e. the propulsion system states affect the airframe 

dynamics and vice versa, indicating that there is adequate airframe/propulsion interaction 

to warrant an integrated control design. Analysis of the control distribution matrix, Bp, 

indicated that the flaps, <5^ and £pvp, and thrust vectoring, &py, are primarily airframe 

controls with very little effect on the propulsion system dynamics, whereas the fuel flow, 

WF, and the area controls, A78 and A8, affect both the airframe and the engine dynamics. 

Thrust vectoring, ^y, is the primary pitch control while the flaps are direct— lift devices 

which provide direct control of the flight path angle. 

For the STOL approach and landing task, it is desirable to provide the pilot with 

decoupled control of flight path and airspeed. Open— loop analyses of the plant indicated 

that it was not possible to provide flight path control independent of the pitch attitude 

control. Therefore, the control law design objective was chosen as that of providing 

independent control of pitch attitude and airspeed from pilot inputs with the flaps used to 

augment the aircraft flight path to pitch attitude (-2) response. The desirable flight path to 

9 

pitch attitude response for the STOL task is 

where (s) represents the Laplace operator. The Level I handling qualities requirements^ 


T(s) « l Jlh 

0 desired s+1 /r^ 
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are that 1 /r^ > 0.7 sec - From the numerator of the open— loop pitch attitude response to 
<^ T y input, the value of 1 /r^ was found to be 0.27 sec -1 indicating that flight path 
augmentation will be necessary to obtain an acceptable system. It was not at all clear as to 
how the LQG/LTR-based DMICS methodology could be used to formulate the problem in 
such a manner as to directly synthesize the desired flight path to attitude response. 
Therefore, prior to applying the DMICS methodology, classical single-input single-output 
control law design techniques were used to augment the flight path response. Since the 
flaps are direct lift devices, a constant gain loop closure from angle of attack to the flaps, 
i.e. ^pL= — and <^p'p= — Kg a with a = C ft Xp, will effectively increase the value of 
1 /tq . The numerical value of matrix C Q is listed in the Appendix. Based on open-loop 
evaluation of the effectiveness of the flaps, the gains and I<g were chosen such that Kg= 
— Kp Based on the theory of coupling numerators^, the root— locus technique was used to 

Q 

study the effect of these loop closures on the zeros of the numerator of (s) transfer 

&TV 

function. A detailed example of this procedure can be found in Ref. [10]. A value of = 3 
deg/deg resulted in 1 /r^= 0.52 sec - *. Although 1 /r^ could be increased further by 
increasing Kp doing so will result in excessive flap deflections. Therefore, K^= 3 deg/deg 
was used as a reasonable compromise between handling quality requirements and flap 
deflection constraints. With the flap loops closed, the control inputs available for the 
feedback compensator design are 

^ = [WF, A78, A8, ^ tv ] T . 

For the rest of this paper, the vehicle model with the flap loops closed as above, will 
be considered to be the design model. 

Feedback Compensator Design and Evaluation 

The design specifications for the feedback compensator were chosen as follows: 

(1) Track airframe and engine commands in a decoupled manner and with zero 
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steady-state error for step commands. 

(2) Maintain "adequate" stability margins in all control loops to guarantee stability in 
the presence of unmodelled dynamics and variation in model parameters. 

(3) Avoid "excessive" control deflections and rates to prevent non-linearities due to 
control deflection and/or rate limiting. 

As stated earlier, the control design objective is to provide the pilot with 
independent control of aircraft airspeed and pitch attitude. In the STOL mode, thrust 
reversing is used to allow higher frequency control of the forward velocity. To effectively 
use thrust reversing, the engine is maintained at Military power level. This power setting 
will use the engine fan speed (N2) and engine pressure ratio (EPR) set points from the 
baseline engine control schedule. Thus the airframe and engine outputs to be controlled are 

y p = [v, q y , N2P, EPR] T 

where 

V = aircraft airspeed (ft/sec) 
q v = pitch variable = q+0.10 

N2P = engine fan speed (% of maximum allowable rpm at operating condition) 

EPR = engine pressure ratio 

The C matrix corresponding to these outputs is listed in the Appendix, and the D matrix 

r i“' 

is 0. Note that the choice of q y above reflects the desire to track pitch rate commands at 
high frequencies and pitch attitude commands at low frequencies. This choice is consistent 
with the work reported in Ref. [4]. 

Based on performance requirements, stability robustness to unmodelled dynamics 
such as actuators and sensors, and the open-loop analyses of control effectiveness, the 
desired control bandwidths were chosen as 1 rads/sec for the V loop, 5 rads/sec for the q y 
loop, 5 rads/sec for the N2 loop and 10 rads/sec for the EPR loop. 

As pointed out in an earlier section, and vividly brought forth by an example study 
in Ref. [11], it is very important to properly scale the system of Eqn. (1) prior to applying 
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the LQG/LTR methodology in order for the methodology to lead to the desired results. For 
the present study, an input/output scaling of the following form was used 



resulting in a scaled system of the form 




x = A J[x + IL u^ ; y^ = C x_ (2) 

P P P Ps Ps ’ J Ps Ps P 

where B = B'S -1 and C = S C . Here, B' is the control distribution matrix 
Ps P u Ps y P P 

corresponding to the control inputs u' , and A' is the plant matrix with the flap loops 

r r 

closed. The numerical values of S and S correspond to normalizing the control inputs by 

u «/ 

their maximum allowable variations from trim and the controlled outputs by the maximum 


values to be commanded, are listed in the Appendix. 

Since the vehicle model has no integrators, and a zero steady— state error for step 
commands is desired, integral control action (see Ref. [12]) is provided by appending four 
integrators to the plant — one in each control channel. The augmented plant has the form 


x = Ax + Bu ; y = Cx (3) 

r P T 1 T 

where x = [Xp, UpJ with Xp and Up § as defined before, and y = yp g . The system matrices 
A, B and C for the scaled, augmented plant are given by : 

\K B n i roi 

A = p Ps ; B = ; C = (C 0] (4) 

0 0 J [ij Ps 

where I is an appropriately dimensioned identity matrix. 

The LQG/LTR procedure is based on solving for the optimal compensator that 
minimizes the performance index 

T 

J = E{ 1 im — /[(z T z + py T Vy) 4- pu T u]dt} (5) 

p T^oo T o 

for a system of the form 

x= Ax + Bu + T? 

y = Cx + vXrj ; z = Hx (6) 

where J and rj are zero— mean Gaussian white-noise processes with identity intensity, y are 
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the measurements available as compensator inputs, and z are the controlled plant outputs. 
V is any positive definite symmetric weighting matrix, and v, p, p and the process noise 
distribution matrix T are the design parameters that are used in the LQG/LTR procedure 
to synthesize a compensator that would meet the desired specifications. The procedure 
consists of first designing the Kalman filter such that the filter loop satisfies the 
performance and stability robustness requirements, and then recovering this loop 
asymptotically by tuning the regulator. The Kalman filter gains Kp and the regulator 
gains I<£ are obtained by solving their respective algebraic Riccati equations^. The 
LQG/LTR compensator has the form 

^LQG/LTR^ = (^) 

Detailed discussion of how to choose the LQG/LTR design parameters is available 

in various application studies of the LQG/LTR methodology, see for example Refs. [13] 

and [14]. For instance, it is shown that if the Kalman filter loop is to meet the design 

specifications, then T and u shall be selected in such a way that the Target Feedback Loop 

transfer function matrix GpQ L (s), defined as 

g fol>> = (-hcfsi-Ar'r (8) 

has the desired loop shapes. From the design specifications stated earlier, the requirement 
on ^"FOL^ * s 

<W b >4^ 1 ’ 5 ’ 5 ’ 10 ] ( 9 ) 

i.e. bandwidth of 1, 5, 5, and 10 rads/sec for the V, q y , N2 and EPR loops respectively, k/s 
behavior at loop cross-over for stability robustness, and large low— frequency gains for 
accurate steady— state tracking of commands. 

For the IFPC control problem being considered here, the controlled variables are 
the same as the measurements, i.e. H = C, and the matrices A, B and C are as defined in 
(4). The singular values of the augmented, scaled plant (Uj[G(ju>)] with 
G(ja;)=C(jud— A) 1 B) are shown in Fig. 2. These indicate that the design plant G(s) has 
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the desired loop shapes, so the choice T = B can be used in (8) to define the Target 
Feedback Loop G>pol( s ) in the synthesis of the Kalman filter gains Kp. The block 
diagram for the LQG/LTR compensator design is shown in Fig. 3. In Fig. 3, G p (s) is the 
transfer function matrix of the vehicle model with the flap loops closed. With p = 1 and V 
= I, in (5), an LQG/LTR compensator that satisfies the desired command tracking 
bandwidth requirements was obtained for the design parameter values v = 0.64 and p = 
10^. Note that although choosing a higher value of p would have led to a better recovery of 
the Kalman filter loop and the associated guaranteed stability margins, doing so would also 
have resulted in larger control bandwidth requirements. The singular values of the loop 
transfer matrix at the plant output (aj[GK(ja>)] with K(ju;) = Klqg/LTR^)) are shown 
in Fig. 4. In Fig. 4, the large loop gains at low frequencies indicate the desirable properties 
of good command tracking and disturbance rejection, and the loop bandwidths are from 1 
to 10 rads/sec providing adequate roll-off at higher frequencies to maintain stability in the 
presence of unmodelled dynamics. Detailed closed— loop performance and stability 
robustness evaluation results with this design compensator are presented in the following. 

The closed— loop system tracking performance was evaluated by studying the Bode 
plots of the closed-loop system frequency response to commanded inputs, and time 
histories of closed— loop response to step commanded inputs. The frequency response of 
velocity V to all the commanded inputs (V , q , N2P and EPR ) is shown in Fig. 5 as an 

^ V £ L L 

example of the Bode plots that were studied. Note that the quantities shown in Fig. 5 are 
scaled quantities (corresponding to y and y of Fig. 3), so that the response magnitudes 
can be compared one-to-one to assess the degree of coupling/decoupling in the closed-loop 
system response. As seen from Fig. 5, the feedback compensator provides accurate tracking 
of velocity commands up to a bandwidth of 1.8 rads/sec (here, bandwidth is defined as the 
frequency at which the response magnitude of the primary commanded variable is —3 dB) 
with insignificant response in velocity to other commanded variables. The bandwidths of 
the q y , N2 and EPR loops were 11, 6.9 and 8.7 rads/sec respectively. The larger than 
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desired bandwidths for tracking the V, q y , and N2 commands are acceptable as long as 
these do not require large bandwidths for control deflections and/or lead to instability in 
the presence of unmodelled dynamics. The discrepancy between the desired bandwidth for 
EPR (10 rads/sec) and the achieved bandwidth (8.7 rads/sec) was not severe enough to 
warrant a redesign. The time-domain analysis of closed-loop system response to step 
commands also showed decoupled command tracking with well-damped response and fast 
rise times for the q y , N2 and EPR responses. Note that apart from good command tracking 
performance, it is also important to make sure that the feedback compensator does not 
result in excessive variations in other variables such as compressor and turbine stall 
margins, high pressure temperature turbine blade etc., which are of interest from the point 
of safe operation of the engine. The responses of these variables were studied for the present 
compensator design, and were found to be within allowable limits for maximum step 
commands in the controlled variables. 

Extensive stability robustness studies, using the singular value approach^’ 
were performed for this feedback compensator design. The effect on closed— loop system 
stability due to including some representative actuator and sensor dynamics was first 
studied. The design plant G(s) was augmented with the following sensor and actuator 
dynamics 

Sensors: N2 = N2 ; EPR = EPR 

s s+20 s s+100 

Actuators: WF = -1!L . _M_ WF„ 

s+10 s+50 c 

A = . - 1 - 00 A for Areas A78 and A8 

s 2 +2(0.5)(30)s+(30) 2 s+100 c 
£ 15 £ 

tv _ “777 ®tv„ 

s+15 c 

to obtain the higher order "truth" model G Q (s). The error between the "truth" model and 
the design model was defined as multiplicative error at the plant output leading to the 
following definitions : 
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G 0 (s) = (I+E(s)]G(s) =» E(s) = [G 0 (s)-G(s)]G *(s) (10) 

The stability of the closed— loop system for the "truth" model is guaranteed if the following 
singular value condition is satisfied for all frequencies w: 

■zII+tGKOu))- 1 ] > *[E(j W )] (11) 

In inequality (11), a refers to the minimum singular value, a refers to the maximum 
singular value, and I+(GK) — ^ is the inverse return difference transfer matrix at the design 
plant output (y). The two quantities in (11) are compared in Fig. 6. From Fig. 6 we note 
that the stability condition (11) is violated in the frequency region 15 to 30 rads/sec. 
However, eigenvalue analysis of the closed— loop system with the "truth" model and the 
design LQG/LTR compensator showed the closed— loop system to be stable, thus pointing 
to the conservativeness of the singular value analysis approach. Clearly, then, there is a 
need to develop improved procedures for stability robustness analysis of multivariable 
systems. 

Lehtomaki et alJ^ have shown that for multi— input multi— output systems, the 
minimum singular value of the return difference matrix, calculated at the plant output or 
input, is a reliable measure of closed— loop system stability robustness to unstructured 
uncertainties occurring at the plant input or output. For the present feedback compensator 
design, the lowest value of the minimum singular value of the return difference matrix at 
the design plant output (o;[I+GK(jw)]) was 0.72 and that at the design plant input 
(a[I+KG(jw)]) was 0.65. These values indicate reasonable stability robustness to 
unstructured uncertainties. (Note that if the loop recovery parameter p was chosen to be 
large enough such that the filter loop is completely recovered, then we will have 
<r[I+GK(ju>)] > 1). However, it is important to realize that these stability robustness 
properties are guaranteed for the loops broken at the points (1) and (2) in the LQG/LTR 
compensator design block diagram of Fig. 3. The actual control system implementation will 
be as shown in the block diagram of Fig. 7, wherein the interface between the physical 
system (G (s) = C (si— A' ) -1 B' for the present case) and the control system (K (s)) is at 

lr r t* t* P 
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points (l 1 ) and (2'). It is at these points in the loop that we need "good" stability margins. 
The lowest value of the minimum singular value of the return difference matrix at points 
(1') and (2'), minimum over io of a[I+K G (jw)] and <r[I+G K ] respectively, were found 

r r ir r 

to be 0.02 and 0.1 respectively. These values indicate very poor stability robustness to 
unstructured uncertainties for the system as implemented in Fig. 7. From this analysis it 
would appear that the claim that the LQG/LTR methodology recovers the guaranteed 
stability robustness properties of the LQR (Linear Quadratic Regulator) problem is 
misleading. The stability robustness properties are recovered , but not at the physical 
interface between the plant and the control system where the effect of the modeling 
uncertainties will actually occur. 

In control system design, a classical measure of stability robustness are the gain and 
phase margins^ Although the above results indicate that the implementation system (as 
in Fig. 7) will have much reduced stability robustness with respect to unstructured 
uncertainties as compared to the design system (as in Fig. 3), it has been shown in Ref. [11] 
that as far as structured uncertainty in the form of loop gain or phase variations is 
concerned, the implementation system will have the same guaranteed multivariable 
stability margins as obtained for the design system. This result is based on a 
nonconservative measure of stability robustness for structured uncertainties, called n or the 
structured singular value, that has recently been developed^’ 1 A detailed mathematical 
discussion of the definition and properties of //, and the form of structured uncertainties 
that can be analyzed using /x can be found in Ref. [18]. For the case of loop gain or phase 
variations, this stability robustness measure is given by 

^[M(jo;)] = mi n £[DM(jcj)D~^] (12) 

DcJ) 

where M is the return difference matrix at the loop point where the gain or phase 
variations are being considered (for example, M = I+KG(ja;) for point (1) of Fig. 3), and J) 
is the set of all real diagonal matrices. In Ref. [19], it is shown that if 
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M M (.M] > /^min ’ 0 < < oo (13) 

for some constant f^ m ^ n < 1, then the nonconservative guaranteed multivariable gain and 
phase margins are given by 

Gain Margin : GM = — 

^min 

Phase Margin: PM = ±2sin -1 (^^^) (14) 

2 

For the special case of gain and phase variations being considered here, the optimal 
* 

matrix D corresponding to the minimum structured singular value /z in Eqn. (12) at each 
frequency u can be obtained using a globally convergent algorithm developed by 
Osborne^. This algorithm was applied to the return difference transfer matrices at points 
(1) and (2) with the design LQG/LTR compensator, and resulted in /* m j n = 0.71 at the 
plant input and /^ m j n = 0.73 at the plant output. Using (14), these results imply that the 
closed— loop system will have guaranteed multivariable gain margins of —4.7 dB to 11.1 dB 
and phase margins of ±42.4 deg for simultaneous gain or phase variations at the plant input 
or the output. These are excellent stability margins and meet the military specifications for 
the operation of flight vehicles. 

The stability robustness of the closed— loop system as implemented in Fig. 7 was 
further evaluated by using the more "classical" approach of "breaking" one loop at a time 
(i.e., one loop open and the other loops closed). An example Bode plot for the loop broken 
at the plant input WF is shown in Fig. 8. As seen from Fig. 8, the closed— loop system will 
remain stable for gain increase of up to 15 dB or additional phase lag of up to 50° in the 
WF loop with all other loops closed at nominal design values. The one— loop— at— a— time 
analysis revealed equally good gain and phase margins for individual variations in the other 
input as well as output loops. 

The one— loop— at— a— time analysis is also useful in studying the integrity of the 
system with respect to actuator and sensor failures, and the actuator bandwidth 
requirements for the control loops. The present design was found to be unstable for failure 
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in the q y and N2 loops, and the A78, A8 and <$ T y loops, thus indicating the need to 
provide adequate redundancy in the corresponding sensors and actuators to ensure fail-safe 
operation. The control bandwidths (defined as the loop crossover frequency, for example u> c 
in Fig. 8) for the WF, A78, A8 and ^ Ty loops were 3.2, 1.7, 2.6 and 6.2 rads/sec 
respectively. These control requirements are very much within the actuator bandwidths for 
the representative actuators considered earlier in the stability robustness study. The 
control deflection and rate requirements for tracking maximum commanded values of the 
controlled variables, obtained from time histories of the control input (Up) responses to 
step commands (y ), were also found to be reasonable and within the maximum allowable 
control deflections and rates. 

Compensator Partitioning and Simplification 

The need to separately implement controllers for the airframe and the engine was 
pointed out in an earlier section. In this section, the mathematical approach to partitioning 
the high— order integrated compensator into separate sub-system controllers is presented 
along with the numerical results for partitioning the LQG/LTR compensator obtained in 
the previous section. The partitioned sub— controllers are further reduced in state 
dimension and simplified for ease of implementation. Some example results are then 
presented to compare the closed— loop system performance and robustness for these 
simplified sub— controllers with that for the full— order LQG/LTR compensator. Note that 
partitioning the integrated compensator and then reducing the order of the sub— controllers 
is intuitively more meaningful than reducing the order of the integrated compensator before 
partitioning. In the latter case, any decoupling inherent in the compensator due to the 
decoupled behavior of the integrated system (plant) might be lost in the process of 
reducing the order of the compensator, thus making the partitioning task more difficult. 
Moreover, reducing the order of the sub— controllers would be easier computationally as 
these would already be of lower order than the integrated compensator. 
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Compensator Partitioning 

The mathematical approach to partitioning the integrated compensator into 
separate (state— decoupled) sub-controllers, presented in the following, was motivated by 
the discussion in Ref. [2] on partitioning the integrated system (plant) into state-decoupled 
sub-systems so as to perform the decentralized, hierarchical control system design. For 
ease of discussion and relevance to the present design study, the procedure is presented for 
partitioning into just two sub-controllers, those for the airframe and the engine. The 
procedure, however, is clearly applicable for the more general case of more than two 
sub-controllers. The direct feedthrough term (D ), for the compensator, has also been 
neglected in the following to keep the equations simple. The extension to the case of the 
non— zero D c matrix is quite straightforward. 

Given an integrated compensator of the form 

\ = A c*c + B c 5 y ; i=C c x c < 15 > 

we want to obtain state-decoupled sub— controllers of the form 


Airframe 

? c. = Vc s + Vy a + Vy, 


Controller ' 

u a = C c. x c. + W ae u e 

(16) 

and 



Engine 

*c, = Vc e + Vy. + Vya 


Controller * 

1 

u e = C c e x ce + W ea u a 

(17) 


— T T — T T T 

with u = [u , u ] and e = [e 1 , e“ ], such that the input-output behavior (u to e 

v V .y a .re \ 


response) of the integrated compensator is matched as best as possible. This desired 


partitioning is shown in the simplified block diagrams of Fig. 9. Note that with such a 
partitioning, the states of one sub— controller do not directly affect the outputs of the other 
sub-controller. Therefore the sub-controllers can be built as separate "black-boxes" and 
integrated in the final implementation with the least number of interconnections. The 
division of the integrated compensator inputs (e ) and outputs (u) into the respective 

•7 
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sub-controller inputs (e and e ) and outputs (u and u ) is determined based on the 
results of the open— loop control effectiveness analyses, conducted at the initialization of the 
control design study. This is coupled with an intimate knowledge of the physical system 
being controlled, i.e. the control inputs u„ into the integrated plant (G(s)) are primarily 
effective in regulating the errors e with secondary response in the errors e , and 

i/ ft i/ c 

similarly for u e - 

The state definition of the integrated compensator has no particular significance in 

terms of physical interpretation. Therefore, without loss of generality, the eqns. (15) can be 

assumed to be in the modal form so that the matrix A is block diagonal. The matrices B 

and C are then measures of the modal controllability and observability respectively for the 

compensator modes. The modal compensator states x c are assigned to x ft or x g based on 

whether they are more controllable by e r or e and/or more observable in u„ or u . This 

ya ye a e 

is accomplished by considering the columns of B and rows of C one at a time and 
comparing the relative magnitudes of the elements in the columns/rows. From eqns. (15) to 
(17), note that all the integrated compensator inputs (e ) are also inputs to the state 

J 

equations of the sub-controllers. This means that the modal controllability for the 
integrated compensator is perfectly matched in the partitioning. Therefore in case of a 
conflict in assigning the integrated compensator states to x„ or x . such as the case where a 

a 6 

modal state is "most" controllable by an element of e whereas it is "most" observable in 

ya 

an element of u 0 , the assignation should be based on observability. After the division and 
an appropriate re-ordering of the compensator modal states, the state-space 
representation for the compensator is as follows 






rB„ B„ lie 
cn C 12 y a 

B B 

C21 C22 e 

L yej 
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( 18 ) 


u 


C„ 

C„ 



a 


Cll 

C12 


Ca 




C 

c 



u e 


C21 

C22 


i X CeJ 


From a direct comparison of eqns (16) and (17) with eqn. (18), we get: ^ Ca = ^ Cl > 
B =B , T =B„ , and A„ =A„ , El =B„ and r =B„ . So the only approximation 
involved in the partitioning is that of approximating the contribution of the states of one 
sub-controller to the outputs of the other sub-controller. Simple algebraic manipulation of 


eqns. (16) to (18) leads to the following solution : 


w = c c# 

ae C 12 c 22 

W = C C# 
ea C 21 c 11 


C c a = C c,r W ae C c 21 

V C c 2 r W ea C c 12 


(19) 

In (19), (*)^ denotes pseudo— inverse of (•). For numerical accuracy, singular value 
decomposition based procedures^ are used to solve for the pseudo— inverse. 


For the LQG/LTR design compensator being considered here, from Fig. 3 the 
compensator inputs are the scaled tracking errors, i.e. e^ = [ey g , e^ e^ 2 p g > e EPR,J 
with e v = V — V , and similarly for other errors with (•) denoting scaled quantities. The 
LQG/LTR compensator outputs are the derivatives of the scaled plant control inputs, i.e. 


u = [WF S , A78 s , A8 S , &j.y ] , and the compensator is of 13 in order, i.e. dim(x c ) = 13. 
Earlier open— loop analyses indicated that Apy is primarily a pitch effector while WF, A78 
and A8 mainly affect the engine control variables (N2 and EPR) and the airspeed (V). 
Therefore, the following sub— controller input/output pairings were chosen for compensator 
partitioning : 


e ya ~ Vs ; \~ 6 TV s 

e ye = [ e y g ? e N 2 P s ’ e EPR g ] ’ u e = A8 S ] (20) 

The compensator modal controllability/observability studies, as discussed earlier, 

resulted in a 5 in order airframe sub— controller (dim(x ) = 5) and an 8 in order engine 

Ca 

sub-controller (dim(x ) = 8). This compensator partitioning was validated by comparing 

c e 
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the LQG/LTR compensator singular values s * n § u ^ ar va ^es 

of the integrated compensator obtained by assembling the partitioned sub— controllers (i.e. 
combining (16) and (17) in the form of (15)), and the singular values of the corresponding 
loop transfer matrices (<r-[GK(jo;)]) and closed-loop response matrices (<t ; [T (ju>)], with 
y(s) = T (s)*y„ (s)). Excellent matching was obtained for all these comparisons, indicating 

S Cg 

that the closed— loop performance and the desirable feedback properties of the integrated 
LQG/LTR compensator are retained in the partitioned sub— controllers. Example 
comparison results are presented after the sub-controller order reduction and simplification 
discussed in the following subsection. 

Sub— Controller Order Reduction and Simplification 

For ease of implementation, it is desirable to further reduce the order of the 
sub-controllers and simplify the structure by eliminating the insignificant feedbacks. Each 
of the sub-controllers, from eqns. (16) and (17), can be represented in the form 


= Vc, + B c&, ; u i = Vci + D ci e ;, 


( 21 ) 


with e^.= col.[e y ., e ], and corresponding definitions of and D V The model reduction 
technique based on internally— balanced realizations^^ can then be used on a system in the 
above form to obtain reduced order sub-controllers. The internally— balanced realization 
approach consists of transforming the system into a set of states such that each state has 
the same controllability/observability index and, furthermore, the states are ordered in 
decreasing value of controllability/observability index. The reduced order system is then 
obtained by eliminating the states that are "less" controllable/observable. Note that in 
general, a frequency— weighted internally— balanced reduction technique has been shown 
to lead to a "better" lower order approximation for controller reduction. However, in the 
present design study, the frequency weighted technique is quite cumbersome to apply as 
the state-space realization of the frequency weighting, suggested in Ref. [23], will itself be 
of very high order. 
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Using the internally— balanced realization approach, the order of the airframe 
controller was reduced from 5 to 2 and that of the engine controller from 8 to 7. The 
system matrices for these reduced order sub-controllers, consistent with the definitions in 
(16) and (17) and corresponding to the input/output pairing of (20), are listed in the 
Appendix. From the data in the Appendix, we note that there is very little direct coupling 
from the output of one sub-controller into the output of other sub-controller, i.e. W and 
W ae are » 0, and also the engine control variable errors (e^p and e ppR ) do not 
significantly drive the states of the airframe sub-controller, i.e. last two columns of 1? 

cl6 

(r (:,2:3)) are « 0. Therefore, the sub— controllers were simplified with W = 0, W = 0, 

d6 a6 6a 

and T (:,2:3) = 0. Note that the fact that such a simplification was possible in the present 
example indicates that the coupling between the airframe dynamics and the engine 
dynamics is not very strong. Such will not be the case in general for the STOVL aircraft 
being considered for the future. 

The loop transfer matrix singular values (<7j[GK(ju>)]) with K(jui) corresponding to 
the compensator obtained by assembling the reduced order and simplified sub— controllers 
are shown in Fig. 10. Comparing Fig. 10 to Fig. 4, we note that the loop shapes obtained 
with the LQG/LTR compensator are preserved in the compensator partitioning and 
simplification. The closed— loop response singular values (^JT (ju>)l) with the LQG/LTR 
compensator and the assembled, simplified sub— controllers are compared in Fig. 11. Note 
that only the maximum and minimum singular values are shown in Fig. 11 to avoid 
cluttering the figure. From Fig. 11, we note that there is no noticeable change in 
closed— loop system performance on partitioning and simplifying the compensators. 

The closed— loop system block diagram with the simplified sub— controllers is shown 
in Fig. 12. The representation in Fig. 12 is consistent with that in Fig. 7, i.e. the 
integrators augmented to the plant for compensator design and the input/output scaling 
matrices are incorporated in the airframe and engine sub-controllers. The blocks marked 
P^, in Fig. 12, correspond to the command shaping filters and are discussed in the 
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following section. Also, as shown in Fig. 12, the engine control variable commands (N2P ge j 
and EPR ge j) will be generated by engine scheduling control system. The discussion of the 
engine scheduling is beyond the scope of the research reported herein, as it requires 
knowledge of the nonlinear behavior of the propulsion system dynamics. 

Pre-Filter Design 

With a satisfactory feedback compensator obtained as above, the next step in the 
DMICS procedure is to design command shaping prefilters to obtain the desired response to 
airframe and engine command inputs. The linear quadratic regulator based prefilter 
synthesis procedure^ was used in the present study. As shown in Fig. 12, separate 
single— input single-output prefilters were designed for each controlled variable as 
decoupled responses to pilot and engine command inputs (£ gt (in.) and V sel- “ d N2P sel 
and EPR ge j respectively) were desired, and the feedback compensator itself provides 
decoupled tracking of the commanded variables. 

The desired controlled variable responses to pilot and engine command inputs are 
listed in Table 1 in transfer function form. In Table 1, the velocity response choice is such 
as to result in a well-damped response with minimal overshoot and a rapid settling time 
with a rise time t r90 = 5 secs, and the pitch rate response (-^-) selection is based on desired 

^st 

short period characteristics for Level I handling qualities^. The N2P and EPR response 

selections are based on the desired command tracking bandwidths. Note that the desired 

q 

pitch variable response ( — ) can be obtained from the desired -3- response using the 

<^st <^St 

definition q v = q+0.10. 

The following steps were followed to obtain the command shaping prefilters that 
will give the desired response for each input /output pair : 

1) The single-input single-output system from the commanded variable to the 
response variable with all the feedback loops closed was first obtained (y.=T.y 
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where is one of V, q y , N2P or EPR). 

2) A reduced order approximation (T r .) of Tj was obtained using internally balanced 
realization. This step is important because the prefilter design procedure results in a 
prefilter of order greater than the "system" for which the prefilter is being designed. 

3) The procedure of Ref. [7] was applied to the reduced order system T to obtain the 
command shaping prefilter P- which will give the desired response. 

4) The order of the prefilter was further reduced using the internally balanced 
realization approach. 

The above process resulted in 4th order prefilters for V, q y and N2P, and a 1st order 
prefilter for EPR. The designed prefilters are listed in Table 2 in transfer function form. 
The Bode plots for the desired pitch rate response and that actually obtained with the 
complete designed control system are compared in Fig. 13. We note good agreement 
between the desired and actual response up to 5 rads /sec — well beyond the short— period 
frequency of 2.24 rads/sec. The other three actual responses, although not shown here, were 
also in agreement with the desired responses for frequencies beyond the specified response 
bandwidths. 

Finally, the frequency response Bode plots of pitch attitude and flight path response 
to pilot stick input, — and respectively, are shown in Fig. 14. From Fig. 14 we note 

^st ^st 

that the flight path to pitch attitude relationship in the region of the short period 
frequency can be approximated by 

6 gp s -1-0. 5 

This is the response that was designed for by considering the angle of attack to flaps ( a -* 
and a -> S^) loop closures prior to the feedback compensator design. 

Conclusions 

The major objective of this study was to gain extensive insight into the Linear 
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Quadratic Gaussian/Loop Transfer Recovery (LQG/LTR) methodology based procedure 
for design of Integrated Flight /Propulsion Control (IFPC) laws. This methodology was one 
of two developed under the Air Force sponsored program on Design Methods for Integrated 
Control Systems (DMICS). Towards this objective, an example of an integrated 
flight/propulsion control system design, mainly based on the aforementioned DMICS 
study, was presented in this paper. The control system design was performed for the 
piloted longitudinal landing task of a modern, statically unstable, fighter aircraft powered 
by a two spool turbofan engine and equipped with a 2D thrust vectoring nozzle. As a result 
of this study, the following conclusions can be drawn regarding the strengths and 
weaknesses of the LQG/LTR based DMICS methodology : 

(1) The major strength of the methodology lies in its ability to address design issues 
such as stability robustness and trade-off between performance and allowable 
control power directly in the synthesis procedure. The procedure for the design of 
the feedback compensator and command shaping is quite straightforward and can be 
accomplished with available control design software. 

(2) The major drawback of the methodology is that it results in a single high-order, 
integrated feedback compensator. Although such an integrated compensator is 
intuitively appealing in that it will lead to "optimum" performance, it is 
unacceptable because of the difficulty of implementation as well as the desire to 
physically separate the flight and propulsion controllers. 

(3) The methodology is mainly geared towards synthesizing command tracking control 
laws and so one of its weaknesses lies in its inability to address the issue of plant 
augmentation. In this example study, classical single— input single-output control 
law design techniques had to be used to synthesize the desired flight path 
augmentation prior to the application of the DMICS methodology. 

(4) The methodology makes extensive use of singular values for synthesis as well as 
analysis. As is well known, and was also shown in this example study, the singular 
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value based robustness analysis tends to be highly conservative. It was also shown 
that because of the scaling (normalization) used in the control law synthesis 
procedure, the LQG/LTR procedure guarantees the stability robustness properties 
at a point in the loop which is internal to the compensator as it will be 
implemented. This point is not the physical interface between the actual plant and 
the implementation compensator where the uncertainties will actually occur. 
Reducing the conservativeness of singular value stability robustness analysis and 
resolving the effect of scaling are issues which warrant further investigation. 

To overcome the major drawback of this DMICS methodology, mainly that of one 
high order integrated feedback compensator (see (2) above), a technique was presented in 
this paper for partitioning the integrated compensator into the desired separate 
sub-controllers for the airframe and the engine. The suggested compensator partitioning 
approach was demonstrated by application to the example study, and led to much 
simplified sub— controllers without any significant loss in performance and stability 
robustness as compared to that obtained with the full order integrated compensator. The 
compensator partitioning approach, however, clearly needs to be validated for a system 
exhibiting higher degree of coupling between the sub— system dynamics than was the case 
in the present study. 

Appendix 


Vehicle System Matrices: 


— 5.893e— 02 1.067e-01 -3.860e+01 -3.184e+01 1.410e-02 3.144e-04 
— 2.659e-01 — 2.665e-01 1.948e+02 -4.599e+00 5.196e-04 -1.578e-05 
— 1.541e-03 7.806e-03 — 1.949e— 01 -4.818e— 04 2.564e-05 9.463e-07 


0 

0 

1.0 

0 

0 

0 

1.427e— 01 

— 9.898e-01 

0 

2.006e+02 

0 

0 

7.782e— 01 

1.542e-01 

0 

0 

— 8.485e— 02 -4.191e+00 

1.518e— 01 

3.008e-02 

0 

0 

— 1.655e— 02 

4.263e— 01 

7.934e-01 

1.572e-01 

0 

0 

— 3.502e— 01 

2.295e— 01 

— 1.005e— 01 

— l.992e-02 

0 

0 

1.096e-02 

3.740e— 02 
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2.599e— 04 3.819e-02 2.251e-03 ' 
— 2.106e— 06 1.826e— 04 — 2.957e— 06 
3.744e— 07 3.668e-05 2.673e-06 
0 0 0 

0 0 0 


6.022e+00 -3.434e+02 1.160e+01 
-5.707e+00 2.716e+01 1.040e+01 
1.155e— 01 — 9.024e+01 8.476e-01 
— 1.036e— 01 — 7.954e+00 -1.068e+00 


B 


P 


3.468e-02 

-4.960e-02 

3.436e-05 -2.055e-01 

6.912e— 02 -4.183e-04 

6.925e-02 

-1.455e-01 

1.234e-08 — 2.936e— 04 

7.104e— 05 

— 5.452e-01 

-8.100e-03 

7.132e— 04 

5.507e-08 1.068e-04 -8.383e-05 

— 7.973e— 02 

0 

0 

0 0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

1.469e-01 0 

0 

0 

0 

0 

5.366e-02 0 

0 

0 

0 

0 

1.813e — 02 — 4.3026+01 

— 2.583e+01 

0 

0 

0 

1.643e-01 0 

0 

0 


C a = [ — 5.888e-02 2.796e-01 0 0 0 0 0 0 0] 



9.797e-01 1.941e-01 0 0 

0 0 5.730e+01 5.730e+00 

0 0 0 0 

0 0 0 0 


0 0 

0 0 

0 8.723e-03 
0 0 


0 0 

0 0 

0 0 

0 6.803e-02 


Scaling Matrices: 

S u = diag [2.0e-04, 0.02, 0.01, 0.1] 


S = diag [0.05, 0.3, 0.2, 3.0] 

J 


Reduced— Order Sub-Controller matrices: 



-5.873e+01 -4.146e+01 
4.146e+01 — 7.861e— 02 



— 2.684e+01 
9.322e — 01 


C„ = [2.684e+01 9.323e-01] 

c a 


r 


ae 


3.209e— 01 2.491e-03 -5.344e-03 
— 1.482e— 02 1.409e-03 -5.143e-03 


W ao = [ 1.113e-04 6.345e-03 -5.262e-03 ] 



— 1.671e — 01 — 1.907e+01 -1.526e-01 
1.737e+01 — 2.023e+01 -1.434e-01 
9.982e— 02 ^.347e-02 -1.213e+01 
-6.165e— 03 -9.159e-02 -8.559e+00 
1.170e+00 — 1.067e+01 -2.406e-01 
5.297e+00 -1.411e+01 -1.630e+00 
8.426e— 01 -5.655e+00 -3.689e-01 


8.728e— 04 -1.040e+00 
— 1.579e— 02 -1.064e+01 
8.596e+00 -1.698e-01 
— 1.856e— 02 -2.981e-02 
— 1.642e — 02 — 1 . 9280- f~0 1 
1.159e— 01 — 1.032e+02 
— 4.437e— 03 -3.562e+01 
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o o o o 



4.518e— 01 — 2.837e— 01 
1.036e+00 7.068e— 02 
6.662e— 03 9.724^-02 
-4.037e— 02 1.264e— 03 


— 1.894e+01 2.991e+01 
-3.599e+02 2.509e+02 
— 2.467e+02 -2.091e+01 


' — 1.948e— 02 -1.872e-02 1.835e+00 ' 


' 2.858e— 03 

2.746e— 01 1.916e+01 -7.368e-01 


— 1.406e— 02 

— 1.128e+01 2.138e— 01 -2.278e-01 


9.420e-01 

-3.723e— 01 5.392e— 02 2.792e-02 

; r = 

1.849e— 01 

8.107e— 02 5.754e+00 -5.444e+00 

ea. 

— 8.075e— 03 

-6.793e— 02 5.501e+00 -2.934e+01 


— 2.003e-02 

— 2.886e— 02 2.508e+00 -4.505e+00 


4.307e— 04 


C 


Ce 


1.812e+00 1.914e+01 1.991e-01 2.549e-02 5.612e+00 
— 1.366e— 01 6.072e— 01 8.751e+00 -3.204e-01 3.561e+00 
— 2.567e— 01 1.056e+00 -7.188e+00 2.705e-01 4.311e+00 


— 1.449e+00 4.761e— 01 
1.908e+01 — 3.284e+00 
2.292e+01 -3.946e+00 



' — 2.408e— 04 ‘ 

; W p = 

— 1.189e— 02 

ea 

9.894e— 03 


Sub— Controller Simplification: 


r ae< : ’ 2:3 > = 0 


W„ = 0 
ae 


W 


ea 
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Table 1. Desired Response Transfer Functions 


Notation : {K(l/r)[C;wJ = K(s+l/r)(s 2 +2Cu; n s+a^)} 

V _ 0.04(3.13) . q ^ 35.12(0.5) 

V [0.89; 0.36] ’ 6 [0.89;2.24] 

N2P _ 5.0 . EPR = 10.0 

N2Psei (5.0) ’ EPR sd (10.0) 


Table 2. Design Prefilter Transfer Functions 
Notation : {K(l/7-)[£;a> n ] = K(s+l/r)(s 2 +2£u; n s+a> 2 )} 

V c _ 0.013 [0.67; 1.37] (5. 20) (830.21) 
v se l [0 . 89 ; 0 . 31 ] [0 . 56 ; 32 . 42] 

q v c _ 3.60(0.10) (0.51) ( 1 .22) (5.61) 

6 st (0)(0.67)[0.89;2.24] 

N2P c = 486.45[0.53;4.05](66.93) . EPR c ^ 1.30(7.67) 
N2P gel (5.0) (13.03) [0.71; 90.50] ’ EPR gel (10.0) 
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FIGURE 4. - LOOP TRANSFER MATRIX SINGULAR VALUES WITH 
LQG/LTR COMPENSATOR <*j [GK L QG/ LT R ( i u) ] ■ 













FIGURE 6. - STABILITY ROBUSTNESS CHECK USING SINGULAR 
VALUES. 
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(a) CENTRALIZED COMPENSATOR. 



FIGURE 9. - BLOCK DIAGRAM FOR COMPENSATOR PARTITIONING. 



FIGURE 10. - LOOP TRANSFER MATRIX SINGULAR VALUES 
(Oj[GK(ju)I) WITH PARTITIONED, REDUCED AND SIM- 
PLIFIED SUB-CONTROLLERS. 


36 





FREQUENCY, RAD/SEC 

FIGURE 11. - CLOSED LOOP SYSTEM RESPONSE COMPARISON 
WITH LQG/LTR COMPENSATOR AND WITH PARTITIONED, RE- 
DUCED AND SIMPLIFIED SUB-CONTROLLERS. 
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FIGURE 12. - DETAILED IMPLEMENTATION BLOCK DIAGRAM WITH SIMPLIFIED SUB-CONTROLLERS. 
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